Intrinsic point defects and volume swelling in ZrSi04 under irradiation 
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The effects of high concentration of point defects in crystalhne ZrSi04 as originated by exposure 
to radiation, have been simulated using first principles density functional calculations. Structural 
relaxation and vibrational studies were performed for a catalogue of intrinsic point defects, with 
different charge states and concentrations. The experimental evidence of a large anisotropic volume 
swelling in natural and artificially irradiated samples is used to select the subset of defects that 
give similar lattice swelling for the concentrations studied, namely interstitials of O and Si, and the 
anti-site Zrsi, Calculated vibrational spectra for the interstitials show additional evidence for the 
presence of high concentrations of some of these defects in irradiated zircon. 

PACS numbers: 61.72.-y,61.80.-x,78.30.-j 
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I. INTRODUCTION 

Large concentration of defects in materials produce 
strains that affect the crystalline lattice structure, which 
can undergo amorphization, change into a different crys- 
talline phase, or remain in the same crystal phase with 
changes in the lattice parameters. In addition to these 
changes, the electronic properties are affected in a way 
that can produce substantial differences in the physical 
and chemical properties of the original material. Radia- 
tion effects in ceramic materials have been subject of a 
considerable number of studies due to their use in nuclear 
fusion applications, and in the stabilization and immobi- 
lization of nuclear high-level wasteiiS. Defect accumu- 
lation and its consequences for the degradation of the 
mechanical properties have to be understood in order to 
control possible leaching, creep and fatigue of these mate- 
rials. A number of issues have to be investigated, such as 
the stabilities of possible defects, their structure, migra- 
tion energies for interstitials and vacancies, displacement 
threshold energies, defect interactions and clustering, and 
their electronic structure. 

Zircon (ZrSi04) has attracted considerable attention 
due to several physical properties of interest, such as high 
permitivity (zirconium silicates are candidates as a gate 
dielectric in future silicon field effect transistorsi^i) , low 
thermal conductivity and high chemical stability. It has 
also been proposed as a ceramic host for the immobiliza- 
tion of nuclear waste due to its high ability to accom- 
modate actinides by substitution of Zr, and its high re- 
sistance to corrosionM. Natural zircon may contain ura- 
nium in concentrations of up to 20,000 ppm, and samples 
with up to 4 billion years have been dated. Thus, even 
if possibly not the best final material for nuclear waste 
encapsulation, it is an ideal system to study the effects 
of radiation damage. 

After a radioactive impurity undergoes an a-decay 
transition, strong damage is produced in the crystalline 
structure. The heavy recoiling atom (tens to hundreds 
of keV) propagates through the material losing energy 
mainly through atomic collisions, thus producing large 
cascade collisions and in consequence, amorphized re- 



gions. The highly energetic a particle ('^5 MeV) that 
propagates in the opposite direction to the recoil atom, 
dissipates most if its energy by ionization processes, and a 
relatively smaller number of atomic displacements mainly 
at the end of its path. 

When the dose of a-decays is not too large there is a co- 
existence between the amorphized (metamict) region and 
the crystalline phase, the latter containing point defects. 
Barriers to recombination cause the accumulation of large 
concentration of these defects, resulting in changes of the 
unit cell volume and shape. This swelling is measured 
from the shifts in the reflection angle in X-ray diffraction 
(XRD) and it is clearly distinct from the one produced by 
the amorphous region, observed by changes in the total 
density. 

Experiments based on natural samples, Pu-doped zir- 
con and ion-beam irradiation, show a unit cell volume ex- 
pansion of up to 5%. This volume swelling is anisotropic, 
giving an increase of ~1.5% along the a5-axis, and ^^2% 
along the c-axis^. Over geological periods the recovery 
of a large fraction of point defects is observed in nat- 
ural zircons (as opposed to Pu-doped samples^), with 
preferential self-annealing of the ab plane^i. 

Characterization of structural changes in the crys- 
talline phase^ show that the dodecadeltahedron ZrOg is 
the less stable unit in zircon, and the defects are expected 
to be more closely related to it, as the Si atoms are well 
confined in the Si04 tetrahedra. Theoretical studies of 
neutral point defects show that the oxygen interstitial 
is the most stable one^. A non negligible concentra- 
tion is expected at thermal equilibrium. The energies re- 
quired to form native defects in zircon were investigated 
by molecular dynamics simulations of primary-knock-on- 
atom collisions with empirical potentials^. It was shown 
that the effective energy required to displace O from its 
crystal position is '^17 eV, whereas it is larger for Si and 
Zr ('^48, and ^60 respectively). 

The previous theoretical studies of point defects in zir- 
con were focused on the energetics of neutral defects. 
However, due to ionization processes produced by the 
a-particle, and the fast that the resulting defect land- 
scape affects the electronic chemical potential, the de- 
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fects can be charged. The accumulation of charged de- 
fects produce electric fields that would change the kinet- 
ics of accumulation and diffusivity. In the present work, 
we use first principles electronic structure calculations to 
study the influence of very high concentration of neutral 
and charged point defects on the crystalline structure 
of ZrSi04. We perform structural relaxations, studying 
electronic properties, as well as some spectroscopic cal- 
culations of a wide variety of defects, forming a quite 
complete catalogue that is compared with experimental 
data. The effects of order/disorder of these defects have 
not been considered in the present study, where periodic 
supercells were used to describe them. 



II. METHODOLOGY 

Calculations of the electronic properties of ZrSi04 
are performed with the self-consistent ab initio siesta 
methodic, using Density Functional Theory (DFT)ii 
within the Local Density Approximation (LDA)^'i. 
The core electrons are replaced by norm-conserving 
pseudopotentialsi^ in the Kleinman-Bylander formic, 
generated with Zr(4s^4p^4d^5s2), Si(3s^3p2), and 
0(2s^,2p'*) atomic valence configurations. The 4p elec- 
trons in Zr were explicitly included in the calculations 
as semicore states, due to the large overlap with the va- 
lence states. The wave functions are described with lin- 
ear combinations of strictly localized numerical atomic 
orbitalsi^. We used a single-C basis set for the semicore 
states of Zr, and double-C plus polarization for the va- 
lence statesi^. The charge density is projected on a real 
space uniform grid with an equivalent plane-wave cut- 
off of 227 Ry, to calculate the exchange-correlation and 
Hartree matrix elements. 

To simulate the defect structures, the host crystal is 
represented by a supercell generated by repetition of 
the conventional tetragonal unit cell (4 formula units). 
The point defect is then introduced inside this supercell 
(adding atoms for interstitial, removing atoms for vacan- 
cies, etc), that have to be of the right size to describe 
the sought large defect concentrations. We use host su- 
percells with 24, 48, and 96 atoms. The supercell with 
48 atoms is a repetition of the tetragonal cell along the 
c axis, and the one with 96, is a repetition along the a 
and b axis. In this way, different concentrations of defects 
can be simulated by changing the number of repetitions 
in the cell. The supercells with 24 and 48 atoms were 
already used by Crocombette in his study of points de- 
fects, and he estimated that unreal physical interaction 
between periodic images of the defects was of the order of 
0.5 eV. For the high concentrations we are interested in, 
these interactions between defects are expected and only 
their periodic arrangement might be unrealistic. The use 
of more than one defect per supercell can alleviate this 
approximation, and will be discussed later. 

Simulating charged defects with periodic boundary 
conditions, requires handling with care the long-range 



coulombic interactions. The unphysical divergence in 
the energy coming from the long range Coulomb inter- 
actions of a periodic array of charges is compensated 
by a uniform electron-charge neutralizing background^i. 
The minimum-energy structure for each defect was ob- 
tained by relaxing with a conjugate-gradient minimiza- 
tion for the forces and stresses. The lattice vectors and 
the atomic positions were allowed to relax until atomic 
forces and stresses from the electronic structure calcula- 
tions were smaller than 30 meV/A and 6 meV/A'^, re- 
spectively. The formation energies of the defects in the 
different charge states are function both of the electron 
chemical potential and of the atomic chemical poten- 
tials for the species involved in the defect, and will be 
discussed elsewhere^-. 

The local vibrational modes (LVM) associated to the 
defects are obtained from the dynamical matrix in the re- 
laxed structure, by finite differences, i.e. computing the 
energies with each atom displaced by ~ ±0.016A along 
each of the cartesian coordinates. The eigenvalues and 
eigenvectors of the dynamical matrix give the vibrational 
frequencies and eigenmodes in the harmonic approxima- 
tion. Previous studies show a deviation of ~5% for vi- 
brational frequencies of point defects in semiconductors. 

The Born-effective charge tensors, Z*^^^, are used 
to determine the infrared intensities, and the splitting 
between longitudinal (LO) and transverse (TO) optic 
modesSS. They are computed by finite differences as the 
change in the i component of the polarization (computed 
with the Berry phase formalism.^ ^) when the atom a is 
displaced along the j direction. 



III. CRYSTAL ZrSi04 

Zircon (ZrSi04) has tetragonal lAi/amd space group. 
The structure consists of alternating Si04 tetrahedra 
and ZrOs triangular dodecadeltahedra sharing edges and 
forming chains parallel to the crystallographic c axis (Fig. 

A body-centred unit cell can be chosen, containing 
four formula units. The structure is fully described with 
4 parameters (see table 0: The lattice constants a and c, 
and the internal parameters y and z. Zirconium and sili- 
con atoms are located in the 4a and 46 Wyckoff positions, 
(0, |, i) and (0, j, |), respectively, and the oxygen atoms 
in the 16h sites (0,y,z). Each O is shared between one Si 
and two Zr atoms, at distances 1.62, 2.12 and 2.25A re- 
spectively. These 4 atoms are in the 0-site's mirror plane 
that can be perpendicular to the a or & axis. The relaxed 
parameters calculated are in good agreement with exper- 
imental values. 

The Born effective-charge tensor is shown in table ^ 
and is in excellent agreement with the results obtained 
with linear response calculations^. Because of the local 
symmetry in the Zr and Si sites, the Born charge tensors 
for these atoms are diagonal, with only two independent 
components, Z| and Z^ parallel and perpendicular to 
the crystallographic c-axis, respectively. The large value 
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FIG. 1: {Color Online) (a) Edge-sharing Si-tetrahedra and 
Zr-dodecadeltahedra along the c-axis, corner sharing in the 
a6-plane. The dodecadeltahedra also have edge-sharing with 
each other along the afe-plane. (b) and (c) show the projection 
of the unit cell in the ab plane, and in a plane perpendicular 
to the 6 axis (the second layer of polyhedra has been removed 
in (c) for a better visualization). The dashed circles in (b) 
show the empty spaces in which we place the interstitials. 



of Z^j, and Z^^ with respect to the nominal ionic elec- 
tronic charge is interpreted in terms of mixed covalent- 
ionic bonding due to hybridization with oxygenSS. For 
O, the Born charge tensor is anisotropic, with the charge 
in the mirror plane being larger than the nominal ionic 
charge (Z=-2), and the tensor component perpendicular 



to the mirror plane being smaller. 

Seven infrared active modes (3A2ii-l-4E„) are predicted 
by group theory in ZrSi04 at F. The induced polariza- 
tion for E„ modes is perpendicular to the c axis, while 
Pl2u modes are observed for electric field vectors paral- 
lel to c. The spectrum between 650 and 1400 cm^^ is 
characterized by Si-0 stretching vibrations (E„(4) and 
A2m(3) modes respectively-'^^) and the 0-Si-O bending 
modes £^(3) and A2n(2). Our computed frequencies for 
these modes are 851, 969, 414, and 602 cm~^ respectively, 
with a maximum deviation of 3.8% with experiments. 
The modes below 400 cm~^ are related to Si04 group 
motions against Zr atoms, and motions of Zr atoms. 



IV. POINT DEFECTS 

We have considered a catalogue of neutral and charged 
intrinsic defects, that includes interstitials (X^) and va- 
cancies {y x) of the three elements {X) present in zircon, 
Zr and Si anti-site defects (Zrs^, and Sizr), Frenkel pairs 
{Xpp)^ and other combinations of interstitials and va- 
cancies. 



A. Volume swelling 

The effect of point defects in the lattice structure after 
relaxation is summarized in Fig. [SJ where we show the 
change in the lattice parameters relative to the neutral 
structure for the defects studied. Some defects produce a 
considerable distortion of the crystalline lattice, whereas 
others such as oxygen interstitial, or the vacancies have 
a much smaller effect. The magnitude of the distortion 
depends on the concentration of defects (supercell size 
in our study). For the values of interest in this work, a 
clear trend is observed, with an almost linear dependence 
between the distortion and the concentration (see Fig. 
|3Jl. For the highest concentration studied (one defect 
per unit cell), we explored the effect of defect disorder 
by introducing two point defects in the supercell made 
from two unit cells. Fig. |21 shows a dispersion in the 
swelling due to the different configurations studied in this 



TABLE I: Calculated structural parameters for crystalline 
ZrSi04 compared to experimental values. 





This work 


Ref. [22] 


Experiment™ 


Volume 


129 (A^) 


127 


131 


a 


6.59 (A) 


6.54 


6.61 


c 


5.96 (A) 


5.92 


5.98 


y 


0.068 


0.0645 


0.0661 


z 


0.184 


0.1945 


0.1953 


d(Si-O) 


1.62 


1.61 


1.62 


d(Zr-0°') 


2.25 


2.24 


2.27 


d(Zr-0'') 


2.12 


2.10 


2.13 


O-Si-O 


96° 


97° 


97° 




116° 


116° 


116° 



TABLE II: Born effective charge tensors for Zr, Si, and O 
atoms. The Born charge for O, correspond to the atom lo- 
cated at (0,«,u), with the mirror plane perpendicular to x. 
For Zr and Si, only the diagonal elements are shown. For 
O, the whole charge tensor is presented, as well as the corre- 
sponding eigenvalues. See also Ref. 22] 



Z*z^ (5.33 5.33 4.66) 

ZJi (3.21 3.21 4.38) 

/-1. 19 \ r-1.19' 

Z5 -3.11 -0.15 -3.16 

V -0.34 -2.26/ [-2.20 
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Defect candidates 

FIG. 2; {Color online) Lattice parameter swelling (in %) for the set of defects considered in this work. An average in the basal 
(ab) plane is done to compare with the relaxed crystalline structure. Xi, Vx,Xfp denote interstitials, vacancies and Frenkel 
pairs of X, respectively. The different charge states are denoted by ±g and the different concentrations of defects by xn, with 
n denoting the number of repetitions of the tetragonal unit cell (cells with 24, 48 and 96 atoms, for n = 1, 2, 4 respectively). 
Unless stated, the neutral defect with n = 2 cell is considered. 



way. Even if the effect is sizeable it does not affect our 
arguments or conclusions, though further work in this 
direction will be needed. 

In general, there is a considerable anisotropy in the re- 
laxation, with different changes in the c-axis and in the 
a6-plane. Experimentally^, an increase in c of ~2% is ob- 
served, while a increases by a smaller amount (~1.5%). 
Based in this evidence, we can focus our study on those 
defects that give a behaviour in the lattice swelling sim- 
ilar to the one observed, eliminating those that produce 
a decrease in the cell parameters, or a larger change in a 



than in c. The defects with the stronger distortions are 
the interstitial of silicon, Si^, and the anti-site of zirco- 
nium in the position of silicon, Zrg^ . The calculations for 
charged systems were done with the supercell of 2 tetrag- 
onal unit cells, repeated along the c-axis. The artificial 
coulomb interaction between cells should be then more 
important in the afe-plane, but we observed that the lat- 
tice parameter c is more strongly affected by the charge 
state than ab. 



B. Structural relaxation and electronic structure. 



TABLE III: Frequencies for IR active modes in ZrSi04 (in 
cm-i). 



Mode 


This work 


Ref. [22] 


Expt-i^ 


A2„(l) 


371 


348 


338 


A2„(2) 


602 


601 


608 


A2„(3) 


969 


979 


989 


E„(l) 


295 


285 


287 


E„(2) 


402 


383 


389 


E„(3) 


414 


422 


430 


E„(4) 


851 


867 


885 



1. Interstitials 

Initially, we place the interstitial atoms in regions with 
empty space (see dashed circles in Fig. and per- 
form conjugate gradient minimization of the forces and 
stresses. We have explored different starting points, 
showing here the relaxed structures with minimum en- 
ergy 

For Zrj and Si^ the interstitial remains along the c- 
axis empty spaces (Fig. lb). Nevertheless, the lattice 
is strongly distorted and the Si04 and ZrOg polyhedra 
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FIG. 3: {Color online) Change in lattice constants (%) as a 
function of the concentration of defects, for interstitials and 
vacancies of the three constituent elements. Results are shown 
for 1 defect for each A'^-repetitions of the tetragonal unit cell, 
with A'^ = 1, 2, 4 and 8. The small arrows show the anisotropic 
swelling in the ab plane for Zri. The different points for 1 
defect/cell refer to different disorder reahzations. 



change to accommodate to the presence of the interstitial. 
There is an anisotropic change in the lattice parameters. 
While Zr; produces large orthorhombic distortion in the 
basal plane (with an increase of up to 3% in the a lattice 
parameter, for a concentration of 1 defect in two unit 
cells, and smaller distortions in the other parameters), 
the Sii distorts the crystal mainly in the c-axis, but also 
in the other directions. Upon removal of electrons from 
Sii, new atomic structures are obtained, with the inter- 
stitial causing oxygen atoms to approach, and neutralize 
the charge around it. In this way, the coordination of 
Si^" increases with n and for Si^"* a six-fold coordinated 
structure results, with Si-0 bonds slightly longer ('--^I.SA) 
than in the initial tetrahedral geometry. In Fig. 21 we 
show the electronic deformation density (difference be- 
tween the solid and atomic electronic clouds), Sp, around 
Sii in two different planes. The geometries for positively 
charged Zr-interstitials are basically unchanged with re- 



spect to the neutral structure. 

In agreement with previous theoretical studies^, we ob- 
serve that the neutral 0^ forms a "dumbbell" -like struc- 
ture (see Fig. |2K), with the axis almost perpendicular 
to the mirror plane of the oxygen in the crystal posi- 
tion. The distances between both oxygen atoms, and the 
neighbouring cations are shown in table IIVI A similar 
structure was also observed in calculations of intersti- 
tial O in ZrOi'^'*. The high ionicity of the crystal forces 
the interstitial to a position where the electronic density 
screens the coulomb interactions. The covalent bond to 
the interstitial oxygen reduces the strength of the Si-0 
bond, and the Si04 tetrahedra loses part of its charge. 
MuUiken population analysis shows that the charges in 
the interstitial and in the lattice oxygen are similar, and 
slightly smaller (6%) than in other oxygen atoms. 

The structure of charged O,^ does not differ substan- 
tially from the neutral defect, with a slight increase of 
the 0~-0 distance (see table HVjl . On the contrary, the 
doubly charged interstitial displaces from the dumbbell 
structure, and forms a bridge between neighbour silicon 
atoms, both having now coordination five. In this struc- 
ture, the relaxed 0,^^-0^ distance (2.19 A) does not ad- 
mit a bond between the oxygens, while Si^— O^^ and 
Si_B-0^"^ do (1.74 and 1.78 A, respectively). This can be 
seen in Fig. [Sji, where Sp around the interstitial oxygen 
differs from the one in the neutral defect (Fig. |SJ:). The 
electronic cloud along the Siy^-OrSis directions (Fig. 
resembles the Si-0 bond in the crystalline configuration. 
The relaxation energy from the initial dumbbell configu- 
ration is of about 1.6 eV, while the energy of the neutral 
defect in the atomic structure of 0~^ is 2.9 eV higher. 
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FIG. 4: {Color online) Atomic configurations for siHcon in- 
terstitial in the -1-4 charged state and corresponding charge 
density-deformation contour plot in the plane defined by Sii, 
Oa and Oc (a), and Sii, Oa and Ob (b). Distances given in 
A and the contour plots are drawn between -0.45 and 0.45e, 
with a step of O.OlSe. Solid (dashed) lines correspond to neg- 
ative (positive) values for the charge. 
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FIG. 5: {Color online) Atomic configurations for oxygen interstitials in the neutral (a) and doubly negatively charged (b) states, 
and corresponding deformation density contour plots in the plane defined by Sia , Oa and Oi (c and d for neutral and charged 
configurations, respectively), (e) shows the contour plot in the plane through SIa, Oi and Sis, in 0~^. Distances given in A 
and the contour plots are drawn between -0.45 and 0.45e, with a step of 0.015e. Solid (dashed) lines correspond to negative 
(positive) values of the deformation density. 



2. Vacancies 

In the Vo , the sihcon atom moves slightly in the direc- 
tion of the missing oxygen, and the remaining Si-0 bond 
lengths are increased by 4%. A new electronic state is 
created in the gap, localized around the vacancy. In the 
charged defect Vq , the Si atom moves back toward its ini- 
tial position, and the Si-0 bonds have deviations smaller 
than 1% with respect to the original Si-0 tetrahedra. 



TABLE IV: Interatomic distances in the dumbbell structure 
(in A), for different concentration of defects (supercells with 
24, 48 and 96 atoms). Charged defect is within the supercell 
containing 48 atoms. 

24 48 96 Ref.f8] O," 



Oi-OA 

Oi-Si 
0,-Zri 
0,-Zr2 
OA-Si 
OA-Zn 
OA-Zra 



1.52 
1.96 
2.29 
2.22 
1.70 
2.24 
2.84 



1.52 
1.69 
2.24 
2.84 
2.00 
2.23 
2.22 



1.52 
1.70 
2.23 
2.82 
1.95 
2.30 
2.22 



1.62 
1.83 
2.24 
2.57 
1.78 
2.20 
2.75 



1.99 
1.69 
2.15 
2.98 
1.80 
2.32 
2.18 



The four oxygen atoms forming the tetrahedra around 
the silicon vacancy remain in their position, with only 
minor deviations towards their centre of mass (~-0.7%). 
For the Vg " configurations, the oxygen tetrahedron ex- 
pands, with the distance from the centre of mass increas- 
ing by 2% for n = 1 and by 5% for n = 2. For the Vzr 
defect a similar effect is observed for the oxygen atoms 
forming the dodecadeltahedra. There is a considerable 
expansion of the Vzr—0'^ bond ('^10%), and a smaller 
change in the length (~-t-4%), that explains the 

anisotropic change in lattice parameters. The expansion 
of the Vzr—O" bond increases for the negatively charged 
vacancies. The Vzr—0^ length expansion decreases (in- 
creases) for negatively (positively) charged vacancies. 



3. Frenkel pairs and anti-sites 

Frenkel pairs of Zr and Si are not favourable because of 
the strong ionic interaction between the interstitial and 
the vacancy. There is a trend for the interstitial to move 
towards the open space of the vacancy, and a restitution 
was observed in the relaxation of several initial structures 
considered. In the case of Opp, the interstitial and the 
vacancy do not interact strongly, and the final configu- 
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FIG. 6; {Color online) Charge redistribution density around 
Zrsi and Si, in the tetrahedral plane. 



ration is similar to an isolated vacancy and a dumbbell 
interstitial, even when both are in the same Si-0 tetra- 
hedra. 

In the Sizr anti-site, the Si atom substitutes a Zr atom 
in the center of the dodecadeltahedron, with a reduction 
of the lengths (-4% for Si-0^ and -6% for Si-O"). In Zrg^, 
the Zr forms a tetrahedron with four oxygen around it, 
with bond lengths of 1.9A(+17% longer than the Si-0 
bonds). The Si-0 distances increase slightly (~-|-l%) 
close to the defect, due to the increase in the lattice pa- 
rameters. Fig. Elshows the contour plot of Sp in the plane 
defined by Zrsi and a pair of oxygens in the tetrahedron. 
In the charged configurations, the bond lengths do not 
change substantially. 



C. Vibrational properties 

The infrared (IR) reflection and absorption spectra 
of zircon is considerably affected by radiation damage. 
Changes in peak frequency, decrease in band intensity, 
and line broadening of the spectra of crystalline zircon 
can be related to changes in the lattice bond lengths. In 
addition, new features appear for highly damaged sam- 
ples, with an increase in the IR intensity with increasing 
dose. One of the new spectral features is the appearance 
of an absorption sharp peak near 796 cm~^ in moder- 
ately damaged zirconS^. This peak is in a gap of the vi- 
brational density of states (VDOS) of crystalline zircon. 
Its origin is unclear, and it has been proposed that it 
could be due to a radiation-induced intermediate phase, 
or to distortions in the boundaries between crystalline 
and amorphized regions^. 

We computed the vibrational properties for intersti- 
tials of O, Si and Zr, and for the anti-site defect Zrsi, in 
the relaxed structures with 48 atoms. IR active local vi- 
brational modes associated to the defects were observed 
for Oi and Si^ in a range of energies corresponding to the 
gap in the VDOS (729 and 743 cm"! respectively). The 
eigenmodes are shown in Fig. [T] They mainly involve the 
stretching of Si-0 bonds close to the interstitials. Other 
localized modes were observed for these point defects, 
but they are not IR active. 



The crystal modes are affected by the presence of the 
defect. The changes are more important in Si^, probably 
due to a stronger interaction with oxygen atoms around 
this interstitial. For 0^, the frequencies of the crystalline 
modes are almost unaffected (though some splittings are 
observed for the low frequency modes) showing that the 
defective structure does not differ much from the crys- 
talline one. In the charged configuration (O^^), the inter- 
stitial is bonded to two silicon atoms, hence considerably 
affecting the vibrational modes of the Si-0 tetrahedra. 
Some weak IR active modes appear in the gap for this 
defect. Zirconium interstitial affects mainly the low en- 
ergy part of the spectra (below 500 cm~^), where the 
Zr-0 interactions are more important. The Si-0 stretch- 
ing modes are softer, probably due the increase in the 
distance between silicon tetrahedra, consequence of the 
large expansion in the a6-plane. For the anti-site defect 
Zrgi, the high energy modes, involving the stretching of 
the tetrahedra bonds, are weaker. The low energy modes, 
mainly related to relative displacements of the Zr and Si- 
tetrahedra sublattices, are less affected. 



V. DISCUSSION 

Experimentally, the change in the lattice parameters 
can be directly monitored using x-ray diffraction. Point 
defects would be randomly distributed through the crys- 
tal lattice. If the number of accumulated a-decays is not 
too large (~ 2 x lO^^a/g), and we assume that ~220 
atoms are displaced by each a particle.^, we can expect 
concentrations of 2 x lO^^defects/cm'^. That is roughly 
one defect for each two unit cells, the kind of concentra- 
tion considered in this work. Similar concentrations have 
been used to study the effect of doping on the lattice pa- 
rameters of semiconductors^, but the changes observed 
are considerably smaller (~0.19%). We performed sim- 
ulations of Si interstitials in the crystalline silicon struc- 





(c) 



FIG. 7: {Color online) Local Vibrational Modes for Oi and 
Sii point defects, (a) and (b) show the eigenmodes for oxygen 
interstitial at 729 cm~^ (IR active) and 870 cm~^ respec- 
tively, (c) shows the IR active mode for silicon interstitial at 
734 cm~^. The relatively high frequency of mode (b) is an 
indication of the bond strength between the oxygens in the 
dumbbell. (Distances in A) 
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ture with a supercell of 64 atoms (roughly the same con- 
centration of defects), giving a volume sweUing of 0.2%, 
much smaller than the values obtained for zircon with 
the same approach. In ceramic materials used for nuclear 
applications radiation swelling is more important [/^V/V 
up to 8% in SiC) and both defect accumulation and ra- 
tio between amorphized and crystalline phases play an 
important role^^. 

When compared with experiments, our calculations at 
those concentrations, single out three candidates as pos- 
sible originators of the lattice swelling: O^, Si^ and Zrgi. 
Some concentration of these three defects should be there 
to account for the experimental findings, even if other 
defects can be present at lower concentrations. In the 
following we discuss the compatibility of this conclusion 
with other facts known about the system. 

For weakly damaged zircon, the fraction of crystalline 
phase is larger than the fraction of amorphized regions, 
and the diffraction maxima are weaker than in undam- 
aged crystals. The anisotropic thermal displacement fac- 
tors have been used to quantify the effect of the radia- 
tion damage in zirconi, indicating that the damage has 
more important effects on Zr and O atoms than on Si, 
which remain inside the Si04 tetrahedral units due to 
the stronger binding of Si with O. In contrast with these 
findings, calculations of threshold displacement energies^ 
(Erf) show that for neutral atoms Erf(0)^E(j(Si)iErf(Zr). 

However, it is important to emphasize that E^ can 
change substantially with the charge state. Charged de- 
fects can migrate through the lattice with a lower mi- 
gration energy than neutral ones. In fact, the atomic 
configurations for oxygen interstitials are similar to the 
ones observed for Zr02 and IIf02^^^- In these materials, 
oxygen is incorporated and diffuses in atomic form, act- 
ing as an electron trap, while neutral defects tend to form 
strong bonds with other lattice oxygens. The negative-U 
behaviour favors the appearance of charged interstitials 
that are less tightly bound to other oxygens. 

In a real material, we can expect all kind of defects 
being present, with the swelling produced by some of 
them being compensated with the deflation due to oth- 
ers. The presence of interstitials requires that vacancies 
are formed somewhere in the lattice, but the effect of 
these vacancies on the lattice parameters is smaller (see 
Fig. 121). The formation of Zrsi leaves a Vzr and a in- 
terstitial of Si. The combination of these three defects is 
compatible with a larger swelling along the c-axis than 
in the a6-plane, as observed experimentally. The forma- 



tion energies for the point defects studied^^, indicate that 
the most probable defects would be interstitials of oxy- 
gen and the anti-site Zrsi. Depending on the chemical 
potentials of each species (/iO: MSi and /izr) interstitials 
of Si and Zr might be present. The IR active mode at 
~730 cm^^ would give further evidence for the presence 
of oxygen, and maybe silicon, interstitials. 

^^Si nuclear magnetic resonance studies of damaged 
zircon indicate changes in the Si local environment to- 
wards a polimerization of the structurc^^ during amor- 
phization. This means that Si— O— Si bonds are formed 
in the amorphous region and the initially isolated Si04 
tetrahedra are connected through an oxygen bridge. This 
polymerization requires that on average the number of 
oxygens per Si is smaller in the amorphized region, 
favouring the presence of interstitial oxygens in the rest 
of the system. As we have seen forms bridging struc- 
tures between Si atoms. But also the presence charged Si 
interstitials can give similar polymerization and change 
the chemical shifts of crystalline Si. 



VI. CONCLUSIONS 

Simulations of high concentration of point defects in 
ZrSi04 have been performed to study the effect of these 
on the lattice swelling under radiation damage. Using dif- 
ferent sizes for the host crystalline supercell, we obtain 
various concentrations of periodically repeated defects. 
A roughly linear dependence between the swelling and 
the defect concentration has been observed for these or- 
dered defective crystal. Based on experimental evidence 
of anisotropic swelling, we have selected a set of defects 
as good candidates to be responsible for the lattice ex- 
pansion in crystalline zircon. Vibrational properties in 
radiation damaged samples are also used to justify our 
conclusions. These include interstitials of oxygen and sil- 
icon, and the anti-site defect Zrg^. 
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